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ABSTRACT 

We report the discovery of twin kHz QPOs in the X ray flux of XTE J1807- 
294, the fourth accreting millisecond pulsar (AMP). This is the second AMP 
exhibiting twin kHz QPOs. In contrast to the first case, SAX J1808. 4-3658, 
the frequency separation Av between the kHz QPOs is consistent with the pulse 
frequency (190.6 Hz), not with half that value, confirming for the first time from 
pulsation measurements the inference, based on burst oscillations, that 'slow 
rotators' (spin frequency less than 400 Hz) have Av approximately equal to the 
spin frequency. While the QPOs move in frequency together over a range of 
more than 200 Hz, Av remains constant with an average value of 205±6 Hz. 
Variability components were found in the 5-130 Hz range similar to those seen 
in other LMXBs. The correlations between the QPO and noise frequencies are 
also similar to those in other sources, but shifted by a factor of 1.59 in kHz 
QPO frequencies, similar to the factor 1.45 shift found for SAX J1808. 4-3658. 
Our results argue in favour of a spin-related formation mechanism for twin kHz 
QPOs and against a spin-related cause of the shift in the frequency correlations. 

Subject headings: binaries: close — pulsars: individual (XTE J1807-294) - 
stars: neutron — X-rays: binaries 



1. Introduction 

The accretion-driven millisecond X-ray pulsar (AMP) family is gradually growing. Since 
the discovery of its first member (SAX J1808. 4-3658, Wijnands & van der Klis 1998) six more 
have been added: XTE J1751-305, XTE J0929-314, XTE J1807-294, XTE J1814-338 and 

1 Astronomical Institute "Anton Pannekoek," University of Amsterdam and Center for High-Energy As- 
trophysics, Kruislaan 403, NL-1098 SJ Amsterdam, Netherlands. 

laboratory for High Energy Astrophysics, NASA Goddard Space Flight Center, Greenbelt, MD 20771. 

3 Department of Astronomy, University of Maryland, College Park, MD 20742. 



- 2 - 



recently IGR J00291+5934 and HETE J1900. 1-2455 (respectively: Markwardt et al. 2002; 
Remillard 2002; Markwardt et al. 2003a; Markwardt & Swank 2003; Galloway et al. 2005; 
Vanderspek et al. 2005). All of these are in transient low-mass X-ray binaries (LMXBs) 
showing faint outbursts (peak count rate 10-100 mCrab) every few years that last a few 
weeks. Their pulse frequency (i.e., the neutron-star spin frequency) lies in the range 185- 
599 Hz and the orbital period is between 40 minutes and 2.5 hours. See Wijnands (2005) 
for an updated observational review. 

XTE J 1807-294 is the fourth-discovered accreting millisecond pulsar. The source was 
detected first on February 13 </l , 2003, in periodic scans of the Galactic bulge region using the 
Proportional Counter Array (PCA) instrument on board the Rossi X-ray Timing Explorer 
(RXTE). These scanning observations show a rising flux, with the epoch of maximum X-ray 
flux (count rate of ~ 58 mCrab) occurring between February 19* h and 22 nd (Markwardt 
et al. 2003b). On February 21 st , 2003, in the first pointed observation of the source by the 
PCA, 190.6 Hz pulsations were detected, confirming this source as an accreting millisecond 
pulsar (Markwardt et al. 2003b). An orbital period of ~40 minutes was determined by 
Markwardt et al. (2003a), still the shortest among AMPs. A pulse analysis of the outburst 
was performed by Markwardt (2004). 

XTE J1807-294 had the longest outburst among the AMPs observed up to now. It 
was followed with the PCA for five months, although after a few weeks the source became 
too weak to measure its aperiodic variability. As noted in Falanga et al. (2005), the X-ray 
flux decay did not show a break to a faster decline, such as in the other AMPs (Wijnands 
2005) and attributed in those objects to the propeller effect (e.g., Gilfanov et al. 1998). 
Falanga et al. (2005) described the energy spectrum as a sum of a disk blackbody and a 
thermal Comptonization component, without evidence of spectral lines or reflection. No 
optical or infrared counterpart has been reported to date. Assuming a distance of 8kpc, 
Campana et al. (2005) placed an upper limit of 4xl0 31 erg on the unabsorbed 0.5-10 
keV quiescent luminosity. 

Studying the power spectra of LMXBs is one of the most direct ways of analyzing 
the accretion flow in these systems. In particular the high-frequency power-spectral features 
constitute a probe of the innermost regions of the accretion disk, where gravitational fields are 
extreme (see van der Klis 2000, 2004, for reviews). The first discovered AMP, SAX J1808.4- 
3658 (Wijnands & van der Klis 1998), in its 2002 outburst provided surprising information 
on the relation between the neutron-star spin and the kHz QPOs and confirmed expectations 
concerning the burst oscillations. While the burst oscillations (cf. Strohmayer et al. 1996) 
occurred at the pulse, and inferred spin, frequency (Chakrabarty et al. 2003), the separation 
Az/ between the kHz QPO frequencies was consistent with half 'the pulse frequency (Wijnands 
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et al. 2003). Although based on only a single simultaneous detection of twin kHz QPOs, 
this result carried important consequences for our understanding of LMXBs, as it falsifies 
simple spin-orbit beat-frequency models, yet requires an intimate link between spin and kHz 
QPOs. This led to several new proposed models for the nature of the interaction between 
disk flow and neutron-star spin (Wijnands et al. 2003; Lamb & Miller 2003; Kluzniak et al. 
2004; Lee et al. 2004) replacing, modifying or adding to earlier proposed models (e.g., Miller 
et al. 1998; Stella & Vietri 1998; Abramowicz et al. 2003). So far, all sources with inferred 
spin frequencies v spin > 400 Hz ('fast rotators') have Au u spin /2 and the ones with 
v S pin < 400 Hz ('slow rotators') have Au ss v apin . 

Several correlations have been found between the frequencies of power-spectral compo- 
nents in Z and atoll sources (the two main types of neutron-star LMXBs; Hasinger & van 
der Klis 1989), as well as among black hole candidate systems. Initially, the WK relation 
(Wijnands & van der Klis 1999) and the PBK relation (Psaltis et al. 1999) were identified. 
Including several more correlated components, this set of correlations was synthesized into 
a 'universal scheme of correlations' among atoll sources by van Straaten et al. (2003). Com- 
paring AMPs to other atoll sources, van Straaten et al. (2005) found that the correlations 
in two of the AMPs are shifted relative to those of the other atoll sources. The simplest 
description of the shifts was that the two kHz QPOs had frequencies a factor ~ 1.5 lower 
than in the other sources. 

In this paper we analyze all the RXTE PCA data on XTE J1807-294, and report the 
discovery of pairs of simultaneous kHz QPOs. The separation Av between the seven pairs 
of kHz QPOs we detect is consistent with the pulse frequency, and confirms for the first 
time from pulsations that the slow rotators have a kHz QPO frequency separation that is 
approximately equal to the neutron star spin frequency. Finally, including also the low- 
frequency variability components in our study, we find a shift in the frequency correlations 
by a factor ~1.5, similar to that seen in SAX J1808. 4-3658. 

2. Observations and Data Analysis 

We used all the RXTE PCA data from XTE J1807-294 available in NASA's High Energy 
Astrophysics Science Archive Research Center taken between February ¥l th and July 29*^, 
2003. 

For the color analysis we used the Standard 2 data, which have a 16-s time resolution 
and 129 energy channels. We divided the channels into four bands in a way that is often used 
for neutron-star sources: band A: 2.0-3.5 keV, band B: 3.5-6.0 keV, band C: 6.0-9.7 keV, 
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band D: 9.7-16.0 keV. The count rates in these precise bands were estimated by interpolating 
linearly between the corresponding energy channels. The background was estimated using 
the standard bright source background models. From the counts in each band, we computed 
the colors and intensity as: hard color=D/C, soft color=B/A, intensity=A+B+C+D. We 
then normalized colors and intensity to the Crab values nearest in time (Kuulkers et al. 
1994) and in the same PCA gain epoch (e.g., van Straaten et al. 2003). 

The timing analysis was done using the GoodXenon data, with the original l/is time 
resolution data rebinned into 1/8192-s bins, and including all 256 energy channels. We per- 
formed fast Fourier transforms (FFTs) of 128-s data segments, fixing the frequency resolu- 
tion and the lowest available frequency to ~0.008 Hz; the highest available Fourier frequency 
(Nyquist frequency) was 4096 Hz. No background subtraction or dead-time correction was 
made prior to the FFTs. The resulting power spectra were "weeded" (the pulsar spike at ~ 
190 Hz was removed at full frequency resolution) and a Poisson level was subtracted. Fol- 
lowing Klein- Wolt (2004) we first estimated the Poisson noise using the Zhang et al. (1995) 
formula and then (after inspecting the ~3000-4000 Hz range and finding no unexpected fea- 
tures) shifted it to match the level between 3072-4096 Hz, where no intrinsic power should 
be present, but only counting statistics noise (this shift was in all cases smaller than 0.15% 
of the previously estimated Poisson level). Then we normalized the power spectra using the 
rms normalization (van der Klis 1995). Several observations with similar power spectra and 
colors were averaged into sets labeled A-H in order to improve the statistics (see Table 1). 
The sets were ordered by increasing kHz QPO frequencies. 

To fit the overall broad-band power spectra we use a fit function consisting of the sum of 
several Lorentzians in the so-called u is max representation", described in Belloni et al. (2002). 
In this representation, if z/ is the Lorentzian's centroid frequency and A its HWHM (half 
width at half maximum) , v max = a/z^ 2 + A 2 gives the characteristic frequency of the feature 
(near the centroid if it is narrow and near the half- width if it is wide). The quality factor 
Q = u /2A is a measure for the coherence of the variability feature. Its strength is given 
by the integral power whose square root, in the normalization we use, is the fractional rms 
amplitude of the variability. Four to five Lorentzian components were needed, all of them 
significant to more than 3<r, with one exception noted in Table 3. In some cases in order to 
avoid a meaningless negative coherence, Q was fixed to zero, which is equivalent to fitting a 
zero-centered (z/ = 0) Lorentzian. From now on, consistent with previous work (e.g., Belloni 
et al. 2002; van Straaten et al. 2003; van der Klis 2004) we refer to these components as 
Li, where the L stands for Lorentzian and V is the label identifying the component ('&' for 
break, V for upper kHz, T for lower kHz, '/i' for hump, 'LF' for low frequency QPO). 
Following this notation we call Lj's characteristic frequency vi and its coherence Q{. In 
order to measure the kHz QPO centroid frequencies and their frequency difference Az/, we 
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performed exactly the same fits to the same data sets, now representing the Lorentzians' 
frequency by their centroid frequency While of course the two representations are 

mathematically equivalent, refitting was necessary to calculate the correct error bars for the 
frequencies. We split the sets in the original observations (see Table 1) and checked wether 
this led to a better measurement of Av. Only in one case the statistical significance of the 
kHz QPOs increased (obs. 80145-01-03-00 in set H; Az/=197.8±7.6 was closer to u spin ) but in 
the final analysis we conservatively used the original datasets. Hereafter, when we talk about 
frequencies this refers always to characteristic frequencies {v max ) as given above; whenever 
centroid frequencies (v ) are intended this is explicitly stated. 

Power law functions were fitted to the z/j vs. v u relations (where V is one of the labels 
mentioned above) in order to quantify the correlations between these frequencies. For this, 
we used the FITEXY routine of Press et al. (1992) that takes into account errors in both 
coordinates and fitted straight lines to the log ^ vs. log v u relations. 

To quantify the shifts in the XTE J1807-294 correlations and compare them with the 
shifts in the other AMPs we used the same method as van Straaten et al. (2005): we took 
as a reference the atoll sources 4U 0614+09, 4U 1608-52, 4U 1728-34 and Aql X-l, and 
included only the frequency sets with v u < 600 Hz (as noted in van Straaten et al. (2005), 
the behavior of the power spectral components above that value becomes more complex; see 
also Section 4). We multiplied the upper kHz QPO frequencies of XTE J1807-294 by a factor 
/, which we varied between 0.01 and 5 with steps of 0.001. In each step we fitted a power law 
to the XTE J1807-294 and atoll source points together. The step with the best x 2 gives us 
the optimal shift factor. We also calculated the shift factors in z/j, following exactly the same 
procedure, but multiplying the V; L instead of u u by a factor /. In the above mentioned atoll 
sources there were only very few detections of L LF . Therefore, to obtain some information 
about the shift in this component, we took as a reference the Llfi component of the low 
luminosity bursters IE 1724-3045 and GS 1826-24 (van Straaten et al. 2005), and again 
applied the same procedure. 



3. Results 

Fig. 1 shows the light curve for all the observations, as well as the evolution in time of 
the soft and hard colors (for definitions see Section 2). As noted in Falanga et al. (2005) 
three peaks in the intensity occurred on March the 9 th , lA th and 25 th ; a peak on March the 
30 t/l was instrumental (caused by a flare in PCU4) and removed in the current analysis. The 
overall trend in the colors is a gradual decline following a slight increase during the first 
two days of observations. The color-color and color-intensity diagrams are shown in Figs. 2 
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and 3, respectively. The region of the color diagrams spanned by the observations, as well 
as the overall shape of the power spectra (see below), is reminiscent of the island state in 
atoll sources (see van der Klis 2004). As can be seen from Figs. 1-3, there is a tendency 
for the colors to decrease with increasing QPO frequencies (from set A to set H), but the 
correspondence is not monotonic. 

Two simultaneous kHz QPOs (see Fig. 4) were detected in seven of the eight datasets 
used in this paper; the best-fit centroid frequencies of these QPOs are listed in Table 2 
together with centroid frequency differences and ratios. As discussed below, the identity of 
the Lorentzian below L u is not clear in set A as it could be fitting either Li or LhHz, or 
both, so we did not include this value in the analysis of the kHz QPO centroid frequencies. 
The measured differences Av between the kHz QPO centroid frequencies are consistent with 
being equal to the pulse, and inferred spin, frequency of 190.6 Hz (see Fig. 5), with no 
evidence for variations over the range the QPOs were detected. The weighted average of the 
seven values for Av was 205.1±6.4 Hz, slightly (2.3a) higher than the pulse frequency. 

The broad-band power spectra of the eight sets of observations can be seen in Fig. 6, 
with their multilorentzian fits. The x 2 /d.o.f. of the fits were between 155/168 and 192/143. 
The best-fit parameters of the power spectral components are presented in Table 3. 

The multilorentzian fits to the power spectra comprise a broad component at low fre- 
quencies (5-10 Hz), L b (see Section 2 for nomenclature); two components at intermediate 
frequencies, one of them narrow at somewhat lower frequency and the other one wider and 
at higher frequency, respectively: L LF (15-40 Hz) and L h (25-130 Hz); and two narrow 
components at high frequencies (the kHz QPOs), Li and L u with characteristic frequencies 
100-370 Hz and 340-570 Hz, respectively. In two of the sets (E and G) Lh was not detected, 
while in set H it has a frequency of ~130 Hz and could also be identified as L^hz ( see van 
Straaten et al. 2002, 2003, 2005, for a description of this component in atoll sources and 
AMPs — the values of the correlation parameters given below do not change significantly 
when calculated excluding in set H). The Li of the two lowest frequency sets (A and B) 
are wider than the rest (Q~0.5 compared with Q~ 1.5) and in set A L e lies slightly below 
the vi vs. v u correlation (see Fig. 7), suggesting that we may be looking at a blend between 
Li and LhHz- I n set A we find a narrow Lorentzian (Q=8) at higher frequency (27 Hz) and 
a wider one (Q=1.4) at lower frequency (16 Hz). We call the higher-frequency one Lh and 
the lower-frequency one Llf] note that this implies Qlf < Qh for this data set. This choice 
has only a small effect on the numbers reported below and does not affect our conclusions. 
In the sets where L h was not detected (E and G), Qlf is also low, which might indicate that 
Lh and Llf are blended. 

The correlations between our measured characteristic frequencies are compared to those 
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in other sources in Fig. 7a. As in SAX J1808. 4-3658, our source shows correlations that 
are shifted with respect to those of the atoll sources. Our results are closer to those of 
SAX J1808. 4-3658 than to those of the other sources, but are somewhat different also from 
that source. The shift factors between the correlations in XTE J1807-294 and those in 
the atoll sources, calculated as explained in Section 2, are presented in Table 4. Note that 
due to the intrindic dispersion present in the frequency-frequency relations the reduced x 2 
(Xr — X 2 /d.o.f.) is larger than 2 in three of the four shifts presented. We take this into 
account when deriving the errors by using A% 2 = Xr- The L b and L h correlations show 
shifts in v u by factors of 1.46±0.05 and 1.72±0.05, respectively. The lower kHz QPO can be 
made to match the atoll sources by shifting it by another factor (~1.08) in u u only, but is 
also consistent with the option that both v u and V£ are shifted by the same average factor 
as derived for Lb and (as van Straaten et al. 2005 reported for SAX J1808. 4-3658). This 
can be seen in Fig. 7b, where v u and v t of SAX J1808.4-3658 and XTE J1807-294 have 
been multiplied respectively by 1.454 (the value found by van Straaten et al. 2005) and 1.59 
(the weighted average of the v\> and factors; see Table 4). The figure clearly shows that 
with these shifts the V£ vs. v u correlation, and for v u < 600 Hz also the i/& and Uh vs. u u 
correlations, match those in the other sources. 

Following van Straaten et al. (2005) and previous works, to prevent crowding of the 
figures, the Llf components, with frequencies usually between Vb and Uh, are not displayed 
in Fig. 7 but their frequencies are plotted vs. u h separately in Fig. 8. The u LF values we 
find in XTE J1807-294 are high, above 10 Hz, like in SAX J1808.4-3658 (van Straaten et al. 
2005) and in 4U 1820-30 (Altamirano et al. 2005), but are somewhat below the power law 
extrapolated from the low luminosity bursters. 



4. Discussion 

We found simultaneous twin kHz QPOs in the power spectra of XTE J1807-294. This 
is the second accreting millisecond pulsar to show twin kHz QPOs. As the spin period of 
the neutron star can be inferred precisely and with high confidence from the pulse frequency 
this provides excellent conditions for testing theoretical models of kHz QPOs (see Section 1; 
van der Klis 2004 for an overview). In SAX J1808. 4-3658 the twin kHz QPO centroids were 
~195 Hz apart, consistent with half the inferred spin frequency [y svin = 401 Hz) and not with 
v S pi n as the original beat-frequency models had predicted (see Section 1). In XTE J1807-294, 
we find instead that the twin kHz QPO separation is near once the inferred spin frequency. 
Au is larger than u spin by 2.3o\ This has previously been seen only in 4U 1636-53 (Jonker 
et al. 2002) and may be difficult to accomodate in beat-frequency (Lamb & Miller 2001) or 
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relativistic resonance (Kluzniak et al. 2004) models. 

In subsequent observations of different sources, a picture is gradually being built up of 
the relations between the burst oscillation frequency u burs t, the pulse frequency u pu i se , the 
kHz QPO separation frequency Au, and the relation of these observed frequencies to the 
neutron-star spin frequency u spin . We currently have two cases where Uburst ~ v vu \ se (and 
none where these frequencies differ), three examples of Au pa Vburst and six examples of 
Au pa fburst/^ (s ee Table 6). However, so far the only direct information about the relation 
between u pu i se and Au was the one case (SAX J1808. 4-3658) where Au pa u pu i se /2. The 
case of XTE J1807-294 constitutes the first example of Au pa u pu \ se . All these findings are 
consistent with the simple picture that (i) the pulse frequency is always the spin frequency, 
(ii) the burst oscillations occur always at the spin frequency, and (iii) for u spin < 400 Hz 
('slow rotators'), Au pa u spin while for u spin > 400 Hz ('fast rotators'), Au pa u spin /2. The 
strict dichotomy expressed in conclusion (iii) may be due to small sample size; some overlap 
between the slow and fast rotator groups might well be found in further observations. The 
model of Lamb & Miller (2003) predicts that kHz QPOs separated by once and half u spin 
can occur simultaneously, but this has not yet been observed. 

In SAX J1808. 4-3658, where twin kHz QPOs were detected only once, several commen- 
surabilities between the centroid frequencies of the kHz QPOs and the spin frequency of the 
form ujfi = \iTLUkfi — nu sp i n \ (where z/^o indicates the centroid frequency of component Lj, j 
and k designate u and I or vice versa and m and n are integers) were consistent with the 
data, with in particular u u $ = \^^e,o ~ ^^spinl a remarkably good match. We searched for 
this in XTE J1807-294, but with seven detections of kHz QPOs moving over a 200 Hz range 
in frequency, no commensurabilities remained consistent with the data as the QPOs moved. 
The data are also not consistent with a constant ratio between kHz QPO centroid frequencies 
(Table 2). This suggests that, because in SAX J1808.4-3658 only a single twin kHz QPO 
pair was detected, it was possible for both the commensurability noted by Wijnands et al. 
(2003) and the 7:4 ratio proposed by Kluzniak et al. (2004) to arise by chance. 

The correlations between the frequencies of the several power spectral components de- 
tected in LMXBs (cf. Section 1) are a relatively new tool to study these systems, and we 
are still in the process of figuring out their significance and even their phenomenology. The 
value that we find for the shift in these correlations (~1.59) is similar to, but significantly 
different from, both the value of van Straaten et al. (2005) for SAX J1808.4-3658 (~1.45) 
and the small-integer ratio 3/2=1.5 (a value that could suggest a relation with relativistic 
resonance models; Abramowicz et al. 2003). Yet, the description of the shifts proposed by 
van Straaten et al. (2005) as the most parsimonious one based on the SAX J1808. 4-3658 
data (all u u and u^ values are too low by the same factor close to 1.5) provides a remarkable 
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match to the data also for XTE J1807-294 (Fig. 7b). In SAX J1808.4-3658 only a single v t 
value was available on which to base the conclusion that V£ was shifted by the same factor as 
v u . In XTE J1807-294, blindly applying the same 1.59 shift (derived from the v b and vs. 
u u correlations) to the seven well-distributed i/£ and v u values, produces a very good match 
to the atoll-source V£ vs. u u relation. 

Altamirano et al. (2005) recently reported a similar shift in the frequency correlations 
of the atoll source 4U 1820-30, but by a factor of ~1.17. This suggests that the mechanism 
responsible for the shifts may also be present in non-pulsing sources and that the shift 
factor in v u may be sometimes far from 1.5. Nevertheless, the small sample size of the 
'shifted' sources makes it too early for final conclusions. Descriptions involving different 
shift factors for different components, and involving shifts of other components than the kHz 
QPOs are still possible and, indeed, formally required by our fit results (see Table 4). A clear 
complication is the nature of the shifted correlation diagram (Fig. 7b) at v u > 600 Hz. There, 
the band-limited noise component labeled L b in XTE J1807-294 follows the L b 2 track of the 
atoll sources, and the Lh track lies above that for the atoll sources or SAX J1808. 4-3658. 
Power law fits to the z/j vs. v u relations of XTE J1807-294 and the atoll sources (Table 5), 
show that the slope of the v b vs. v u relation in XTE J1807-294 is much less steep than in 
the other atoll sources, and actually consistent with that of the u b2 points (see Fig. 9), and 
that the v h points define a steeper relation in XTE J1807-294 than in the other sources. 
The overall impression created by all of this in Fig. 7b is that of a bifurcation in both the 
v b and the u h vs. v u relations when v u exceeds 600 Hz. Perhaps such a bifurcation could 
be related to the increasing separation between general-relativistic azimuthal and radial 
epicyclic frequencies when approaching the compact object (e.g., Abramowicz et al. 2003), 
or to the onset of some form of rotational splitting (cf. Titarchuk & Muslimov 1997), but at 
this stage this must remain conjectural. 

Further analysis of other AMPs and atoll sources is needed to gain a better insight 
into the frequency correlation shifts. Only a few useful power spectra were obtained up to 
now for XTE J1751-305, XTE J0929-314 and XTE J1814-338 and no twin kHz QPOs were 
seen in them. The recent second outburst of XTE J1751-305 (Grebenev et al. 2005; Swank 
et al. 2005) or forthcoming outbursts of other new or known AMPs can help to give us a 
more complete view of this shift phenomenon. In any case, our results for XTE J1807-294 
when combined with those for SAX J1808. 4-3658, strongly support a spin-related mechanism 
for kHz QPO production and strongly suggest a spin-independent cause of the shift in the 
frequency correlations. The former is due to the first detection of twin kHz QPOs separated 
by the spin frequency in a confirmed slow rotator, and the latter to the fact that similar 
shifts are observed in the frequency correlations of at least two AMPs with very different 
spin frequencies. 
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Table 1: Observations used for the timing analysis. 



Set 


ObsID 


Start date 


Det. a 


PDS b 


Count rate (s v ) c 
Total - Per PCU 


Bkg. (s- 1 ) d 
Total - Per PCU 


A 


80145-01-02-00 


2003Mar05 


3-4 


126 


266 - 70 


102 - 27 


8D1 A K D1 D/1 HQ 
oU140-Ul-U4-Uo 


zuuoiviarro 




80145-01-02-03 


2003Mar06 










B 


80145-01-02-02 


2003Mar06 


3-5 


99 


308 - 78 


108 - 27 


80145-01-02-01 


2003Mar06 




rdi 4^ m rn 














sni m ni 
ouiio-ui-ui-uo 


zuuoiviai u i 










C 


80145-01-01-04 


2003Mar01 


4-5 


104 


391 - 91 


115 - 27 




80145-01-01-00 


2003Mar02 










D 


80145-01-01-01 


2003Feb28 


5-5 


67 


399 - 80 


135 - 27 


80145-01-03-02 


2003Mar07 


E 


80145-01-02-06 


2003Marll 


2-4 


124 


215 - 74 


84 - 29 




70134-09-02-00 


2003Feb27 










F 


70134-09-02-01 


2003Feb27 


2-5 


169 


267 - 84 


89 - 28 




80145-01-02-04 


2003MarlO 










G 


80145-01-03-01 


2003Marl2 


3-4 


87 


238 - 70 


92 - 27 


H 


80145-01-03-00 


2003Mar09 


2-4 


154 


261 - 77 


92 - 27 


80145-01-02-05 


2003Marl3 



a Minimum and maximum number of active detectors during the observations. 

6 Total number of averaged power density spectra (PDS), each one made from a data 
segment of 128 seconds. 

c Time averaged count rate not corrected for the background. 

d Time averaged background count rate estimated with the standard model for bright 
sources, version 2.1e, provided by the PCA team at GSFC. 
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Table 2: Twin kHz QPQ centroid frequencies with their difference and ratio. 



Set 


v u ,o (Hz) 


vt,o (Hz) 


Av a (Hz) 




B 


352.8±3.5 


127.1±28.0 


225.7±28.3 


2.8±0.6 


C 


374.0±1.9 


184.0±10.5 


190.0±10.7 


2.03±0.12 


D 


393.4±3.1 


195.0±14.3 


198.4±14.6 


2.02±0.15 


E 


449.6±7.1 


216.5±23.9 


233.1±25.0 


2.1±0.2 


F 


466.2±3.0 


246.3±19.0 


220.0±19.2 


1.89±0.15 


G 


490.6±4.6 


253.9±22.83 


236.7±23.3 


1.93±0.17 


H 


563.6±4.6 


359.9±15.8 


203.7±16.5 


1.57±0.07 
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Table 3: Multilorentzian fits to the power spectra. 





Parameter 




Llf 


L h 


L e 


L u 


x 2 

d n f 

U.U.I. 




V ma x{Hz) 


5.3±0.6 


15.6±0.6 


26.8±0.5 


106±30 


337±10 


180 
148 


A 


Q 

rrn Q ( Or, \ 
i may /(J J 


(fixed) 


I. 4±0.3 

II. 6±1.2 


5 4+1 


0.4±0.3 
15±2 


2.8±0.6 
i a 7+1 k 

J_ KJ. 1 _l — L . O 




Vmax{Hz) 


5.3±0.7 


14.45±0.17 


24.7±1.6 


163±23 


354±4 


169 
165 


B 


Q 

rrtW% 

I lllo^ / (J J 


0.26±0.08 
1 1 c_i_r) q 


y -3 
4 Q+1 


0.73±0.18 
i c q_i_i k 


0.6±0.3 
16±2 


5.5±1.0 
15.8±1.2 




Vmax(Hz) 


6.7±0.8 


17.4±0.4 


40±3 


191±9 


375±2 


192 
143 


C a 


Q 


(fixed) 


2.0±0.5 


0.5±0.2 


1.7±0.6 


6.3±0.5 




1 IIIO^ / V J 




9. Oil. 3 


16.2±1.9 


12. Oil. 4 


17.3±0.5 




V ma x{Hz) 


5.7±0.7 


20.2±0.4 


44±12 


202±11 


395±3 


129 
121 


D 


Q 


0.45±0.16 


3.6±1.2 


(fixed) 




5.1±0.6 




rms(%) 


8.8±1.4 


7.8±1.2 


20.3±0.9 


12.1±1.9 


18.9±0.8 






7.5±1.4 


26.1±1.6 




238±25 


449±9 


132 
111 


E 


Q 


0.2 (fixed) b 


1.0±0.5 




1 2+ 1 ' 2 


7 Q+4-4 




rms(%) 


14.2±1.7 


15±3 




17±3 


15.2±1.8 






9.4±1.4 


27.6±0.5 


69±23 


259±16 


465±2 


202 
175 


F 


Q 


(fixed) 


3.2±1.0 


(fixed) 


1.7±0.7 


6.7±0.6 




rms(%) 


13.6±1.5 


8.6±1.2 


20.3±1.1 


13±2 


18.7±0.7 




V m ax{Hz) 


7.8±1.1 


32.1±1.5 




273±19 


492±5 


155 
168 


G 


Q 


0.36±0.11 


1.0±0.2 




1.3±0.6 


8±2 




rms(%) 


14.2±1.2 


16.9±1.4 




20±2 


17.3±1.6 






10.2±0.8 


40.3±1.1 


129±19 


370±18 


565±5 


151 
142 


H 


Q 

rms(%) 


0.22±0.06 
16.0±0.6 


1.7±0.3 
13.7±1.1 


1 2+ 1 - 4 

nt 4 3 c 


2.0±0.8 
17±3 


6.6±1.4 
17.1±1.4 



a An extra component could be added to the fit at ~0.1 Hz, without significant change in the other 
Lorentzians. 

6 Q fixed to 0.2 for the stability of the fit. 

This component was 2.2<7 single trial but when added to the fit function it gave a 2.9cr improvement of 
the chi squared (according to the F-test for additional terms), without significantly changing the rest of the 
components. 



-14- 



Table 4: Shifts in the frequency-frequency relations in XTE J1807-294. 



Component 


X 2 /d.o.f. 


z/j factor 


v u factor 


T o. 


22 / 17 


0.303±0.031 


1.725±0.109 


u 


68 / 22 


0.335±0.030 


1.462±0.048 6 




44 / 20 


0.249±0.017 


1.715±0.048 6 




185 / 21 


0.895±0.083 


1.082±0.069 



"This shift was done using the Llfi of the low luminosity bursters as a reference (see 
text and Figure 8) 

6 The weighted average of these two values is 1.59±0.03. 

C A double shift in this relation (u u and v t equally shifted) gave a factor 1.31±0.26, with 
a x 2 /d-o.f. of 185 / 21. 
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Table 5: Power law fits to the frequency-frequency relations. 



Ui = N x v u 


a 


X 2 /d.o.f. 


Ui — Nx v u a 


a 


X 2 /d.o.f. 


XTE J1807- 


-294 




Atoll sources 






Vb 


1.29±0.23 


4/6 




3.072±0.093 


15 / 14 


VLF 


2.169±0.070 


14/6 


Vb2 


1.53±0.39 


9/11 


Vh 


3.34±0.33 


9/4 




2.530±0.074 


29 / 14 


Vl 


1.64±0.14 


4/6 




1.402±0.026 


181 / 14 



a Five values measured by van Straaten et al. (2002, 2003) coincide with the Lf, track 
at u u > 600 Hz; these were not included in the measurement of the slope of vs. u u . 
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Table 6: Approximate relations between the kHz QPO separation Au, Vburst and v pu i se in the 
sources where at least two of these frequencies have been measured. See Table 2.5 in van 
der Klis 2004 for a more extensive version and complete references. 



Source 


"burst (Hz) 


"pulse (Hz) 




Al// Impulse 


References 


Slow rotators 


XTE J 1807-294 




191 




-1 


This paper. 


4U 1915-05 


272 




-1 




Galloway et al. (2001) 


XTE J1814-338 


314 


314 






Strohmaycr ct al. (2003) 


4U 1702-43 


330 




-1 




Markwardt et al. (1999) 


4U 1728-34 


363 




-1 




Mendez et al. (1999) 


Fast rotators 


SAX J1808.3-3658 


401 


401 


-1/2 


-1/2 


Wijnands et al. (2003) 


KS 1731-260 


524 




-1/2 




Wijnands & van der Klis (1997) 


Aql X-l 


549 




-1/2 




van der Klis (2000)° 


4U 1636-53 


582 




-1/2 




Wijnands et al. (1997) 


SAX J1750.8-2900 


601 




-1/2 




Kaaret et al. (2002) a 


4U 1608-52 


619 




-1/2 




Mendez ct al. (1998a,b) 



"Marginal detection upper kHz QPO. 
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G 
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Time (MJD) 



Fig. 1. — Intensity, hard color and soft color (defined in text) versus time for all the obser- 
vations. Symbols identify the observations used for the timing analysis as indicated. The 
peak count rate (~ 58 mCrab) occurred near February 21 s * (MJD 52691). No strong timing 
features were detected after March 23 th (MJD 52721) because the count rate was too low. 
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T 




Soft color (Crab) 

Fig. 2. — Color-color diagram with a zoom into the area used for timing. Symbols identify 
the observations used for the timing analysis as indicated. The arrow points at the first 
observation and the solid line indicates the path followed in time. 
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Fig. 3. — Color intensity diagrams. Symbols identify the observations used for the timing 
analysis as indicated. Typical error bars are shown. 
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500 1000 1500 

Frequency (Hz) 



Fig. 4.— Twin kHz QPOs in the power specrum of obs. 80145-01-03-00 (in dataset H). Note 
the pulsar spike just below 200 Hz. 
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170 

350 400 450 500 550 
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Fig. 5. — kHz QPO frequency separation Av plotted vs. upper kHz QPO centroid frequency 
v u0 . The thick line indicates the pulse frequency, the thin line the average separation and 
the letters above the frame the datasets used. 
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Fig. 6. — Power spectra in power x frequency vs. frequency representation together with the 
respective fit functions and their Lorentzian components. The pulsar spike was removed 
before rebinning in frequency. 



-23- 




0.1 1 1 1 1 — 1 ■ ■ ■ — — 

100 v u (Hz) 1000 100 v u (Hz) 1000 

Fig. 7. — Characteristic frequencies of the power spectral components (except vlf, which 
is analyzed separately) plotted versus the characteristic frequency of the upper kHz QPO. 
The black symbols correspond to XTE J1807-294, the open ones to SAX J1808.4-3658 and 
the grey ones to the atoll sources . Symbols identify the different components as indicated, 
and the letters above the frame the approximate location of XTE J1807-294 data. The 
observed frequencies are plotted in (a) while in (b) the XTE J1807-294 points were shifted 
by a factor 1.59 and the SAX J1808. 4-3658 points by a factor 1.454 in both v u and vg. Note 
the apparent bifurcations in the L b and L h correlations in (b) for v u > 600 Hz. 
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Fig. 8. — Frequencies of Llf versus the frequencies of Lh, for the sources indicated in the 
plot. The solid line shows the power law fit to the low luminosity bursters IE 1724-3045 and 
GS 1826-24. The dotted line represents a power law with the same index and a normalization 
factor half of that of the solid line. 
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Fig. 9. — Characteristic frequencies of L b and L b2 for the atoll sources 4U 0614+09, 4U 
1608-52, 4U 1728-34 and Aql X-l and for XTE J1807-294, plotted versus v u . The fits to 
the relations are shown as a solid black line (XTE J1807-294), a dashed black line {y#i of 
the atoll sources) and a dashed grey line (f& of the atoll sources). Five points with doubtful 
identification were excluded (see Table 5). 
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